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many works have followed implementing such a count- 
ing H, IE Q with finite cut-offs or proposing a counting 



I. INTRODUCTION 



Effective field theories (EFT) are a powerful tool to 
deal with non-perturbative low energy physics. Over the 
last years, they have provided promising results as re- 
gards a systematic and model independent understand- 
ing of hadronic and nuclear physics. The scale separation 
between long and short distance physics makes the de- 
velopment of a systematic power counting possible. Af- 
ter the original proposal of Weinberg's to design a 
power counting based on applying ChPT to the potential 
ay works 
lii 

in the renormalized S-matrix [5|, |^ which has also been 
pursued to NNLO 't'I . The relation of both the Wein- 
berg (W) and Kaplan-Savage- Wise (KSW) counting has 
been understood as perturbative expansions about in- 
frared fixed points in the limit of small and large scatter- 
ing lengths iSi] respectively (see also Ref. 9] for a discus- 
sion on long range forces in that context). For systems 
with a large scattering length as it turns out to be the 
case in low energy NN scattering, the Weinberg count- 
ing may be modified to iterate the scattering length to 
all orders, but then the connection to ChPT must be 
given up On the other hand the KSW counting, 

although systematic, does not converge at NNLO 7]. In 
Ref. (11| a new counting (BBSvK) involving also the chi- 
ral limit should be invoked. According to these authors, 
one should treat non-perturbatively the NN potential in 
the chiral limit and consider finite pion mass corrections 
perturbatively on top of that. For a recent and more com- 
plete review on these and related issues see e.g. Ref. [T^ 
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and references therein. 

In order to properly define a truly EFT three essential 
requirements must be met. Firstly, one needs a power 
counting scheme, i.e. a dimensionless expansion parame- 
ter, which controls a priori the accuracy of a calculation 
and provides an error estimate of the neglected terms. 
The second requirement is the mathematical need for 
a regularization method and subsequent renormalization 
scheme which is consistent with the physical power count- 
ing, independent on any ultraviolet cut-off. Finally, there 
is the question of practical convergence, which can only 
be decided a posteriori on the light of practical calcula- 
tions. 

The issue of regularization and renormalization in the 
present context is not at all trivial, particularly if the 
power-counting scheme involves summing up some infi- 
nite set of diagrams. There are well known examples 
in the literature that shows that not all regularization 
methods comply to the physical power counting, like e.g. 
■kN scattering for relativistic baryons in dimensional reg- 
ularization 13J. Such an approach is not consistent with 
the non-relativistic limit in the case of heavy baryons. 
One has to either device a specific scheme, the so-called 
infrared regularization j or to make a non-relativistic 
limit first and introduce dimensional regularization after- 
wards jl^. In the context of NN scattering it has also 
been shown |l6j that the renormalized scattering ampli- 
tude for a theory without pions with a contact and a 
derivative term, i.e. a truncated potential, depends on 
whether one uses a cut-off regularization or dimensional 
regularization as an intermediate step of the calculation. 
The problem was latter understood by using other sub- 
traction schemes than the MS in dimensional regular- 
ization, the so-called power divergence scheme (PDS) Q 
(see also Refs. 0, for off-shell (OS) subtraction 
schemes.) 

As usual in field theory, one essential ingredient in the 
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whole construction is the choice of an appropriate regu- 
lator, which may eventually be removed. In perturbation 
theory, there are many such regulators, like e.g. dimen- 
sional regularization. Beyond perturbation theory, like in 
the two body scattering problem one can use those regu- 
lators within an order by order analysis aiming at finite 
renormalized scattering equation. For potentials which 
are purely short range, i.e., contact terms and deriva- 
tive corrections there-off, this problem has successfully 
been tackled. The problem becomes subtle when, in 
addition, long range forces are added. Typically, these 
long range potentials develop some strong singularity at 
short distances. The paradigm of the problem is proba- 
bly best exemplified by the ^So channel with One Pion 
Exchange (OPE) potential. On the one hand, the reg- 
ular OPE alone produces a finite scattering amplitude. 
On the other hand the short range potential is infinite 
but renormalizablc, and hence a finite scattering ampli- 
tude can be built according to the standard principles 
of renormalization theory. This suggests, as it turns out 
to be the case, that the problem may still be renormal- 
izablc within a scheme where the long range potential 
is treated perturbatively. Whether or not this mathe- 
matical requirement is physically justified has been the 
subject of much debate in the recent past, with a negative 
answer after the findings of Ref. [rauJl; iterated OPE 
contributions to scattering observables are not small in 
the triplet '^Si Di channel. The reason may be found 
in the singular and attractive 1/r^ nature of the tensor 
contribution to the OPE potential. 

The mathematical problem actually arises when both 
short and long range potentials are added and none 
of them can be considered small; there seems no ob- 
vious way to renormalize the scattering equation non- 
perturbatively. By this we mean regularizing the equa- 
tions first and removing the regulators afterwards, or at 
least make the mass scale of the regulator much larger 
than any other mass scale in the problem. The main 
interest in making such a non-perturbative renormaliza- 
tion is that any perturbative scheme can be thought of as 
an approximation to renormalized equations. The non- 
perturbative renormalization of NN interaction in the 
singlet ^So and triplet ^S'l — ^£>i channels has been stud- 
ied several times in the literature by different regular- 
ization methods. In Ref. ^ subtraction method was 
developed for the Lippmann-Schwinger equation to con- 
struct a finite T matrix for contact interactions added 
to OPE. Renormalization is indeed achieved by taking 
the subtraction scale to be much larger than any other 
mass scale and checking for independence of results in 
this limit. The resulting description of the ^Sq phase shift 
is only valid to very low energies, requiring for inclusion 
of derivative terms. Unfortunately, the method has not 
been extended to that case. Derivative interactions can 
be included within a tree-dimensional cut-off regulariza- 
tion 20, 21 1 . Actually, in the cut-off regularization analy- 
sis of Ref. i21ij a strong breaking of angular momentum is 
observed unless either a vanishing bare mixing coupling is 



chosen or a (according to the authors) unlikely fine tun- 
ning sets in for a non-trivial fixed point. We will show 
below that indeed the bare mixing vanishes but not in an 
uniform way, but rather following a limit cycle pattern. 
Inspired by the N/D method, the work of Ref. 22] makes 
a re-summation of the KSW amplitudes introducing a on- 
shell potential to which chiral counting is applied. A new 
parameter, which is considered to be of order zero in the 
chiral counting is introduced. The dynamical origin of 
this parameter is unclear. The non-perturbative dimen- 
sional regularization of the OPE potential with derivative 
interactions has been studied in coordinate space as 
well as in momentum space |24j in the ^5*0 channel. In 
both cases a three-parameter fit can be achieved with no 
explicit two pion exchange contribution up to CM mo- 
mentum as large as p '--^ 400MeV. Ref. 24] opens up 
some hope as how the triplet '^S'l Di channel might 
also be renormalized non-perturbatively within dimen- 
sional regularization, but so far there are no practical 
calculations. The early finite sharp cut-off momentum 
space treatments of Ref. Q have been improved by im- 
plementing a better regularization scheme j25j which al- 
lows to make take larger cut-off values, yet finite. More 
recent works supporting the W-counting have also ap- 
peared [2II27I]. 

From a diagrammatic point of view momentum space 
treatments based on the Lippmann-Schwinger equation 
are more natural within a Lagrangian framework and al- 
low explicit consideration of nonlocal potentials. On the 
other hand, the long range NN potentials making use 
of chiral symmetry constraints are local, and for those 
the analysis of non-perturbative renormalization in co- 
ordinate space becomes much simpler, as will be shown 
along this work. In addition, the Schrodinger equation is 
a second order operator and mixed boundary conditions 
define a complete and unique solution of the scattering 
problem in the whole space at both sides of the bound- 
ary. This sharp boundary separation of the space is natu- 
rally formulated in coordinate space for a local potential. 
Boundary conditions for NN scattering were used many 
years ago (see e.g. Ref. and references therein), and 
there has been renewed interest motivated by the devel- 
opments within EFT [UlsHaJ. Actually, the thorough 
analysis of Ref. Isij shows that in the absence of long 
range forces a low momentum expansion of the poten- 
tial within EFT framework for the Lippmann-Schwinger 
equation is completely equivalent to an effective range 
expansion (ERE) and also to an energy expansion of a 
generic boundary condition at the origin in coordinate 
space for the Schrodinger equation. Moreover, the refer- 
ence partial wave analysis of the Nijmegen group 
uses this method to successfully describe the a large 
NN scattering data base, when long range potentials are 
used. While in the first works phenomenological poten- 
tials where used, more recent studies consider potentials 
deduced from ChPT theory with a rather satisfactory de- 
scription of the experimental scattering data j33[ . The 
minimal boundary radius which can still provide an ac- 
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ceptable ^jDOF is about Rs = 1.4 - 1.8fm. Obvi- 
ously, if the radius cannot be lowered without spoiling 
the quality of the fit, the short distance cut-off becomes 
an indispensable parameter of the theory, which cannot 
be removed. The corresponding momentum space cut-off 
A = 2tt /Rs ^ 600MeV is comparable to the one needed 
in early momentum space treatments 0. Within the 
spirit of an EFT it would actually be more appropriate 
to take instead larger A's or equivalently shorter i?s's 
and to check for insensitivity of results in the low energy 
regime. Thus, there arises the natural question whether 
in fact this EFT procedure can be implemented. 

In our previous work 34] we showed that for the 
singlet channel with OPE the boundary radius can be ef- 
fectively removed without spoiling a good description of 
the corresponding phase shift up to the a priori expected 
CM momentum of A: ~ where the Two Pion Exchange 
(TPE) effects should start playing a role. The first or- 
der differential equation satisfied by the boundary con- 
dition of the problem defined in the interval R < r < oo 
as a function of the boundary radius was very helpful, 
since the whole problem could be mapped into a variable 
phase equation |35( of a truncated potential in the region 
< r < R with a non-trivial initial condition at the ori- 
gin, encoding the short distance physics. In this way, the 
long range pions could be eliminated and the evolution of 
the threshold parameters as a function of the boundary 
radius could be determined non-perturbatively. Actually, 
a trivial ultraviolet fixed point limit for the scattering 
length was found non-perturbatively. Remarkably, this 
behaviour coincides with the one found in Ref. [lol | within 
a perturbative treatment. This trivial fixed point at the 
origin implies a fine tuning of the short distance physics 
in order to reproduce the physical scattering length. In 
this paper we want to extend our results for the interest- 
ing case of the triplet "^Si —^Di channel. The solution of 
the boundary condition problem requires solving a cou- 
pled set of Schrodinger equations. Instead of doing so, 
we prefer to directly compute the change of the bound ary 
condition by an equivalent variable phase approach | 35| 
with non-trivial initial conditions which encode the short 
distance physics [s^]. This provides, in addition, a di- 
rect and quite transparent connection to renormalization 
group ideas [EQ- 

Related works in spirit to the present are those of 
Refs. 0,113 and 0. In Refs. 0,113 a square well po- 
tential is used to regulate the short distance behaviour 
simulating a smeared delta function. The renormaliza- 
tion group flow for the potential strength is not uniquely 
defined. This phenomenon is also found in the theory of 
self-adjoint extensions of the Schrodinger operators |37| . 
In Ref. a delta shell regulator located at a finite dis- 
tance is assumed as the short distance potential whereas 
the long distance piece is solved exactly using a distorted 
wave basis. This formalism has been so far used to the 
study of renormalization of repulsive singular potentials 
(like l/r^). The reason may have to do with the need 
for a well defined renormalization at the origin. A com- 



mon feature of both regularization schemes is that the 
wave function at the origin is uniquely determined by 
the regularity condition, u{0) = 0. The boundary condi- 
tion regularization that we use in this paper provides a 
uniquely defined renormalization group flow |34| , to treat 
both repulsive and attractive singular potentials [13 . In 
addition, the boundary condition admits a simple physi- 
cal interpretation: it can be mapped into a variable phase 
shift problem '^ith a truncated potential. This inter- 
pretation directly provides the non-perturbative renor- 
malization flow of low energy parameters and a quite 
transparent analysis of both infrared as well as ultraviolet 
fixed points and limit cycles |38j . 

In this paper we analyze precisely how the energy de- 
pendent boundary condition must change as we move the 
boundary radius for fixed energy to achieve independence 
of physical observables such as scattering phase shifts. 
By doing so we are effectively changing the Hilbert space 
since the wave function in the outer region is defined only 
from the boundary to infinity. An advantage of this pro- 
cedure is that we never need to invoke off-shellness explic- 
itly; at any step we are dealing with an on-shell problem. 
In addition, we work directly with finite quantities and 
no divergences appear at any step of the calculation when 
the boundary radius is taken to zero from above. 

Another advantage of our construction, as it will be- 
come clear along the paper, is that we only need the 
potentials and physical threshold parameters as input of 
the calculation (the cut-off dependence is removed com- 
pletely, so this is not a parameter). This implies, in par- 
ticular, that given this information we never have to make 
a fit (except perhaps for the determination of the thresh- 
old parameters); our calculations are predictions for the 
phase shifts that are consistent, by construction, with a 
low energy expansion up to a given order. Thus, the po- 
tential danger of compromising the low energy fit due to 
a global fit up to SOOMeV may be precluded from the 
start. This is at difference with the standard way of pro- 
ceeding where the low energy parameters are fitted to 
the phase shifts and the threshold parameters are then 
recomputed. Actually, our analysis is equivalent to mak- 
ing a fit only in the low energy region, where explicit 
pions do not contribute, and predicting the intermediate 
energy region. We believe this is a possible and practical 
way of learning about the role of explicit pions in the NN 
interaction. Actually, our motivation was partly to see 
whether OPE can actually be seen in low partial waves 
in the intermediate energy region < k < m^. 

In the present paper we analyze the OPE potential 
{U = 2ijy and ji = AIn/2 ) which reads 
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(4) 



Where Mjv is the nucleon mass, ryijr the pion mass, f.^ the 
pion weak decay constant and gA the nucleon axial cou- 
pling constant. In the numerical calculations below we 
take Mn = 938.92 MeV, = 93MeV , = 138 MeV 
and gA = 1-25. Note that the singularity at the origin of 
the tensor potential 



Ut 



MlNgl 



(5) 



is independent on the pion mass m-^. 

The plan of the paper is as follows. In Sect. |n] we 
present the basic object of our analysis, the variable S- 
matrix which we supplement with general mixed bound- 
ary conditions in the general case of coupled channel scat- 
tering. We also discuss the role played by the irregular 
solutions for singular potentials in the spirit of an ef- 
fective field theory. After that we rewrite in Sect. IIIII 
the variable S-matrix equation for the variable iiT-matrix 
and i?-matrix in a way that the low energy limit may be 
taken. As a result, we find the boundary radius evolution 
of threshold parameters. As a first application we apply 
in Sect. llVl the obtained equations to determine the low 
energy threshold parameters from well established NN 
potentials. In Sect. we study the short distance be- 
haviour of the threshold parameters. There we show that 
one has for the ^5*0 and ^Si —^Di channels an UV fixed 
point and a UV limit cycle for the scattering lengths. In 
Sect. IVll we present our numerical results, both for the 
threshold parameters as well as for the ^Sq and ^Si^^Di 
phase shifts. Finally in Sect. IVlH we present some final 
remarks, conclusions and perspectives for future work. 



II. VARIABLE S-MATRIX WITH BOUNDARY 
CONDITIONS 

In order to generalize to triplet states the results of 
Ref. 34] for the singlet channel case, we introduce the 
variable S-matrix formalism for the general coupled chan- 
nel case. For potentials which are either regular or singu- 
lar repulsive at the origin the procedure is standard |33 
and it has many variants. For completeness and to make 
the exposition more self contained we present here our 
particular derivation which also applies to singular at- 
tractive potentials and at the same time introduce our 
basic notation for the rest of the paper. Although in the 
case under study we are interested in, at most, two cou- 
pled channels, the formalism can be developed for the 
general case with almost no additional effort. 

The scattering amplitude for NN scattering can be 
written as a partial wave expansion 



f = 



1 

2ik 



^(S-l)PKcos0), 



(6) 



where S is the S-matrix for coupled channels. The cou- 
pled channel Schrodinger equation for the relative motion 
reads 



u"(r) + 



U(r) 



121 



u(r) — fc^u(r) , 



(7) 



where U(r) is the coupled channel matrix potential, P = 
diag(/i(/i -I- 1), ... , iNilN + 1)) is the angular momentum, 
u(r) is the reduced matrix wave function and k the CM. 
momentum. We assume for u(r) the mixed boundary 
condition ^ 



u'{R)+Lk{R)u{R)=0, 



(8) 



where Lfc (i?) is a real hermitean matrix in coupled chan- 
nel space, which in our framework encodes the unknown 
physics at distances r below the boundary radius R. In 
addition, we assume the asymptotic normalization con- 
dition 



u(r) Ui„(r) - Uout(?-)S , 



(9) 



with S the standard coupled channel S-matrix. The cor- 
responding out-going and in-going free spherical waves 
are given by 

UoutW = diag(/i+(fcr),...,/i+ (fcr)), (10) 
Ui„(r) = diag(/i-(fcr),...,/i- (fcr)), (11) 

with h^{x) the reduced Hankel functions of order 
hf(x) — xH^^i^{x) ( = e^" ), and satisfy the free 
Schrodinger's equation for a free particle, 

-"outW + ^Uout(0 = fc2uout(0: (12) 



'Uin(^) + —Vi\n{r) = fc^Ui„(r) , 



(13) 



The boundary condition, Eq. (jS)), for the outer boundary 
values problem, Eq. ((TJ and Eq. can be interpreted in 
simple physical terms of a complementary inner problem 
where the potential U(r) acts in the interval R < r < oo. 
If we switch off the potential above a given boundary 
radius R we have, at the boundary 



Lk{R) = u'{R)u-\R) 

X [Ui„(i?)-Uout(i?)S(i?)]"' 



(14) 



where S(R) is the S-matrix associated to the potential 
U(r) acting in the region < r < R, which inherits the 
dependence on the chosen boundary radius R. The equa- 
tion satisfied by the variable S-matrix can be obtained 



1=0 



^ This is the most general boundary condition that makes the cou- 
pled channel Hamiltonian self-adjoint in the interval R < r < oo. 
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from Schrodinger's equation applied to the matrix L(i?) 
yielding 



LkiR)' + Lfe(i?)2 = U(i?) + ;^ - fc' ■ 



(15) 



From here ^ it is straightforward to obtain the equation 
for the variable S-matrix, 



2i.^ 
dR 



S(i?)hW(ii')-h(-)(i?)l U(i?) 



h(-)(ii')-h(+n-R)S(i?) 



(16) 



This is a first order non-linear matrix differential equa- 
tion which can be solved by standard means, provided 
the S-matrix is known at one given scale. One of the in- 
teresting aspects of this equation is that there is no need 
to invoke any off-shellness; for any value of the boundary 
radius we have a different on-shell scattering problem. 
In appendix we show an alternative derivation based 
on continuous deformations of the potential with a fixed 
boundary condition. As we will discuss below, Eq. H16|l 
describes the renormalization group flow of the S-matrix 
as a function of the distance scale R where the long range 
potential is cut-off. 

In the case of a regular potential, Eq. (|16|l has to 
be supplemented with an initial condition at the ori- 
gin, namely the trivial one (corresponding to the absence 
of a potential), and its asymptotic value yields the full 
5— matrix ; 



S(0) 



S = S(c») (regular) (17) 



In this paper we are concerned with the OPE potential 
which has a singular 1/r behaviour at the origin in the 
^5*0 singlet channel and singular 1/r^ behaviour at the 
origin due to the tensor force in the "^Si Di triplet 
channel. While in the single channel the singularity is a 
mild one in the sense that there still exists a unique reg- 
ular solution at the origin, u{0) = (like in the Coulomb 
potential) , in the triplet channel both linearly indepen- 
dent solutions to Schrodinger's equation vanish at the 
origin, and the regularity condition u(0) = does not 
uniquely specify the solution. 

The point of view we take in the present work is that 
of an EFT; low energy physics should not depend on the 
detailed knowledge of the interaction at short distances. 
This applies, in particular, to the case of a singular po- 
tential as will also become clear below. Following the 
lines already sketched in our previous work [34| , we take 
instead the value at infinity as the initial value for the 
variable S'— matrix. Of course, for a short range poten- 
tial, this procedure corresponds to start integrating at 



^ An alternative derivation of Eq. 1151 , closer in spirit to the 
renormalization group and the Callan-Symanzik equation will 
be presented elsewhere Ref. l38l . 



sufficiently large distances (where the potential may be 
neglected). An advantage of this procedure is that by 
construction a unique solution S(i?) is obtained. Even 
for a regular potential, it is clear that a generic choice of 
S(oo) cannot yield by integrating towards the origin the 
result S(0) = 1 besides the very exceptional cases which 
accidentally correspond to the regular solution at the ori- 
gin. Thus, we expect in general an admixture of both 
the regular and irregular solutions, which corresponds to 
a mixed boundary condition close to the origin. 



lim {u'(i?) 



L{R)u{R)} = 



(18) 



In the case of a singular potential both solutions vanish 
and we equally have a unique mixed boundary condition 
as in Eq. (|18|l . Thus, we may define the short distance 
S-matrix as the extrapolation to the origin of a given 
solution at infinity. 



S, = lim S(i?) 



S = S{oo) (general) (19) 



Actually, the precise meaning of the previous limit will 
be the main topic of the present work. We anticipate 
already that we will find ultraviolet fixed points for the 
singlet ^Sq channel and limit cycles for the "^Si Di 
triplet channel. Eqs. (|16|l and (|20|l are well known in 
potential scattering ( for a review see e.g. Ref. js^), but 
they have always been used assuming the trivial initial 
conditions S(0) = 1. 

Obviously, if one would literally use the full S— ma- 
trix and integrate downward, nothing could be achieved, 
since that would correspond to eliminating the full po- 
tential. A more interesting perspective, already pursued 
in Ref. [sj for the singlet channel, consists of re- 
garding the low energy limit of the previous equations, 
extracting the threshold parameters at short distances 
by integrating downward from their experimental values 
and integrate back upward the variable 5'-matrix equa- 
tion to infinity. Physically, this procedure corresponds to 
explicitly separate the OPE contributions on top of any 
good low energy approximation, like e.g. the effective 
range expansion. 

In the case of one channel, like the ^5*0, the S-matrix 
can be parameterized as Si{k,R) = exp(2iSi{k, R)) with 
(5;(fcj_i?) the variable phase. Eq. 1)1611 becomes rather sim- 
ple [35] for s— waves, yielding 

dSoik^ R) ^ 1 ^ 
dR k 

and the obvious conditions both at the origin and at in- 
finity must be satisfied 

lim 5o{k, R) = 5^{k) lim (5o(fc, R) = Sa{k) . (21) 

R—*-0 R.^oo 

The OPE potential in the coupled '^Si—^Di in the triplet 
channel space is given by 



U(r) 



f Us{r) Usd{r) 
Usd{r) Ud{r) 



(22) 
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where 



Us = Uc Usd = 2V2UT Ud = Uc-2UT. (23) 

The two coiiplccl channels S-matrix can be represented 
in the Blatt-Biedenharn (BB or Eigen phase) parameter- 
ization 



S = 



cose 



■ sine 



sin e cos e 







cos e sm e 
- sin e cos e 



(24) 



which will be used along this paper. The relation to the 
standard coupled channel K-matrix is given by 



S = {K + ik) (K-iA:)"^ 



where 



tan(2e) 

- tan 6- 

- tan 6+ 



2K 



12 



i^ll - K22 
1 



^^11 + K22 + 

Ku + K22 - 



K 



11 



K' 



22 



cos 2e 
K11-K22 



cos 2e 



(25) 

(26) 

(27) 
(28) 



Due to unitarity of the S-matrix in the low energy limit, 
fc — > we have 



(S-l);,_; =-2ia;,,ifc''+'+i + ., 



(29) 



with ain the (hermitean) scattering length matrix. The 
low energy limit acquires its simplest form in the Stapp- 
Ypsilantis-Metropolis (SYM or Nuclear bar) parameter- 
ization 



le 



s^i^i cos2e 
i(5i+«2)sin2e 



igi(<5i+'52) sin2e 
e^'^^ cos2e 



(30) 

(31) 
(32) 



which is related to the BB phase shifts by 

5i+S2 = S+ + S^, 
■fx X \ tan(2e) 

""(-^^-^^^ = w^- 

The low energy limit in the SYM representation becomes 

Si —a^k , 62 —a2k^ , e —ao2k'^ . (33) 

The scaled K-matrix, K, has a good low energy behaviour 
and is defined by making an energy dependent transfor- 
mation 



K = fcDKD . 



(34) 



with D = diag(fc'% . . . , fc'-^'). The scaled K-matrix ad- 
mits the coupled channel analog of the effective range 
expansion 

1 



K 



-t- -rk^ + vk^ + 



(35) 



where a, r and v arc the scattering length matrix, effec- 
tive range and curvature parameters respectively. 



III. EVOLUTION OF LOW ENERGY 
PARAMETERS 

In order to take this low energy limit and corrections 
there-off, we introduce the variable or running if— matrix 

S(i?) = (K(i?)-hifc)(K(i?)-ifc)"\ (36) 
as well as the reduced Bessel functions 

ji (x) = xji (x) , yi (x) = xyi (x) , (37) 
i.e. jo{x) = sinx , yo{x) = —cosx. Thus, 

-y = 1 (h(+) + h(-)) . (39) 

Then, we get 

K'{k,R) = (^{k,R%kR)-Y(kR)^\5{R) 

X (^^j{kR)K{k,R)-y{kR)Y (40) 

The scaled K-matrix, K(i?), has a better low energy be- 
haviour and is defined by making an energy dependent 
transformation 



K{R) = fcDK(i?)D , 
with D = diag(fc'i , . . . , fc'" ). We get 



(41) 



K'{k,R) = (^K(i?,fc)-j(fci?)D-i -y(A;i?)Dj U(i?) 

X Qj(/ei?)D-iK(ii,fc)-y(fcii)D^ . (42) 

The scaled K-matrix admits the analog of the effective 
range expansion 

K{R) = -a(ii)-^ + ^r{R)k^ + v(i?)fc^ + ... . (43) 

where a(i?), r(i?) and v(i?) are the corresponding run- 
ning scattering length matrix, effective range and curva- 
ture parameters respectively. In this form the low energy 
limit can be easily taken. Defining the matrix functions 
and their low energy expansion 

AkiR) = ^^D-i = Ao + fc^As + fc*^4 + . . . , 



Bk {R) = y{kR)-D = Bq + ^^^2 + fc*B4 + • • • , 



(44) 



we get the system of coupled equations ^ 
A[a(i?)]-i = -([a(i?)]-iAo + Bo)U(i?) 



^ The Equation for V2 in the coupled channel case is too long to be 
reproduced here, but can be obtained in a straightforward way. 
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dR 



X (Ao[a(i?)]-i+Bo) 
r(i?) = ([a(i?)]-iAo + Bo)U(i?) 



(45) 



X (r(ii)Ao + 2[a(ii')]-iA2 + 2B2) 

+ (r(i?)Ao + 2[a(i?,)]-iA2 + 2B2)U(i?) 

X ([a(i?)]-iAo + Bo) . (46) 

These equations generalize to the coupled channel case 
those already found in Ref. '33| and have to be supple- 
mented with some initial conditions, at e.g. infinity, 



a(oo) = a, r(oo) = r 



(47) 



For the case of s— wave one channel scattering Eq. 142|l 
becomes 



dK{k,R) 
dR 



= UiR) 



, sin kR 

K(k,R) — hcosfci? 

k 



(48) 



where 

K{k,R) ^ k cot S{k,R) , (49) 
yielding at low energies an effective range expansion, 

kcotS{k,R)^ \— + l-ro{R)k^+v2{R)k^--- , (50) 

ao[R) 2 

where 



^ = U{R) (ao - Rf 



(51) 



dR 

dV2 

'dR 



2U{R)R^ 
U{R) ' 



R 



2 1 



R 

"0 



1 - 
R 



R_ 



Sao 

^ro 
12 R 



ro 
R 

- 1 



R 

3an 



- 1 



(52) 
(53) 



V2_ 
i?3 



1 R 
120^ 



1 

24 



These equations have been studied by us in Ref. 
analyzing the OPE in the singlet ^Sq channel. 



(54) 



for 



In the ^Si Di coupled channel case the threshold 
parameters matrices are 



ao 
vq 

V02 



aQ2 
a2 

ro2 
r2 

V2 



(55) 
(56) 
(57) 



The explicit form of the equations for the "^Si —^Di run- 
ning scattering lengths reads 



15R^"' 



02 



225i?'^a!; 



9Udal2 + (ao - R)R^ [("o - R)Us + 6ao2C/.d] 
-15ao2i?'' (-ao + R)Us+ R^ (45ao2 - (ao 

2250^2 R^Us - 30ao2i?' (-45a2 + i?^) Usd ^ 



r 



R) {~45a2 + R^)) Usd 
(~45a2 + R^'f Ud . 



3ao2 ("45a2 + R^) Ud , 



(58) 



Note that all three running low energy parameters ao, 
ao2 and a2 (the explicit R-dependence has been sup- 
pressed for simplicity) are coupled due to the mixing po- 
tential Usd- Thus, it would be inconsistent to take any 
of them as a constant; exact renormalization group in- 
variance requires mixing between the S and D channels. 
As we see the mixing is related both to a non-vanishing 
of the mixing potential Usd and a non vanishing value of 
asd at a given point. If by some accident both vanish at 
a given point, the mixing will vanish. 

The evolution of the low energy parameters can be 
translated into the corresponding evolution of the short 
distance boundary condition as a function of the bound- 
ary radius. Defining the dimensionless quantity 



Cfc(i?) 



RUiR) 



Ru',{R)uk{Ry 



(59) 



and using Eq. ((T3|l we get 

RC'kiR) = Ck(l - Ck) + U(i?)i?2 + F - k^R^ . (60) 

Expanding into powers of the momentum k one gets 

CkiR)^Co{R) + k^R^C2iR) + ... (61) 

For the singlet ^5*0 channel we have, in particular, the 
following relation 



Co 



aoiR) 
R-aoiR) 



(62) 



Note that for i? — s- oo we have a fixed point behaviour 
Co ^ unless a = oo in which case Co — > 1. The evo- 
lution of the boundary condition with the short distance 
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boundary radius for the ^Si—^Di in terms of the running 
scattering lengths is given by 



Triplet ^5';^ -3 Di with mixing Nijmll (Reid93) 



C° = 1 



R{R^ - 45a22) 



45a22) 



45ag2 



= 3 



5(fl-aoo)i?5 
45a^2 + ("00 



(63) 



(64) 



(65) 



i?)(i?5 -45a22) 

(R-dependence has been suppressed for simplicity). 
Again, for R oo we have for non exceptional values 
of theparameters C° ^ , (7°^ ^ and C)^ ^ -2. In 
Ref. [Sg a more detailed study on these issues will be 
carried out. 



IV. DETERMINATION OF LOW ENERGY 
PARAMETERS AND THE THEORY WITHOUT 
EXPLICIT PIONS 

An essential ingredient of our formalism is to parame- 
terize the scattering data directly in terms of low energy 
threshold parameters, like a, r and v, defined through 
Eq. Unfortunately, besides a and tq in the singlet 

and triplet channels, the PWA data base 's^ does not 
provide values for them. They could be obtained from 
a fit to the NN data base in the pertinent channels, at 
sufhciently low energies. Such a procedure turns out to 
be numerically unstable, particularly for the v parame- 
ter, because it depends very strongly on the energy win- 
dow chosen for the fit (see Appendix IB|| . On the other 
hand, the NN data base provides explicit potentials, some 
of them local like the Nijmll and Reid93 potentials, for 
which the variable phase approach may directly be ap- 
plied. In such a way we can uniquely and accurately 
determine all the needed low energy threshold param- 
eters by integrating Eqs. H46|l upwards from the origin 
to infinity with trivial boundary conditions. For illus- 
tration purposes the evolution for the Nijmll potential 
for the threshold parameters is depicted in Fig. We 
remind that these curves represent the low energy thresh- 
old parameters corresponding to a potential truncated at 
a given distance, R. The behaviour at the origin has to 
do with the strong repulsive core of the potential in the 
s— wave channel. Also, the divergence of the scattering 
lengths at about R ^ 2fm signals that a bound stated 
has appeared. The values at infinity correspond to the 
physical values. Our results, in appropriate powers of 
fm, can be summarized as follows for the Nijmll and the 
Reid93 (in brackets) potential 

• Singlet ^5^0 Nijmll (Reid93) 

ao = -23.74(3), ro = 2.67(75) , W2 = -0.48(9) (66) 

• Triplet ^Si without mixing Nijmll (Reid93) 

ao = 5.001(3) , ro = 1.833 , V2 = 0.131(41) (67) 



5.419(22) 1.647(6) \ 
6.504(453) ) 

/ 1.833 0.404(12) 
^ ^ \ - -3.522(66) 

/ -0.131(41) -0.274(64) \ 
~ ' - -3.70(80) J 



(68) 
(69) 

(70) 

(71) 



The ^Si channel without mixing parameters have been 
obtained from the '^S'l Di channel for the ^Si compo- 
nent, and ofQ = l/(a~-'-)oo, complying to the low energy 
expansion of the scaled K matrix, Eq. (|35|l . Although 
we will be using the Nijmll parameters, we have also 
presented the ones corresponding to the Reid93 case to 
provide an idea on the size of errors. 

Once the threshold parameters have been determined 
we can use the coupled channel effective range expansion, 
Ea. (|35|l to find out to what extent does this expansion 
apply. On theoretical grounds we expect this expansion 
to converge within the region of analyticity of the S- 
matrix, which presents a left cut at fc = ±im^/2. In 
Fig. lI2Jl we compare the quality of the ERE including 
LO, NLO and NNLO contributions to the original data 
of Ref. |32| ■ As we see, to describe the data within the 
ERE approach up to the convergence radius m^/2 one 
has to go at least to NLO. The description of the data 
below TO7r/2 is improved, as expected, with higher orders 
in the ERE. Above this region, where OPE should play 
a role, this is not necessarily so. Actually, we see that 
in the ^5*1 end Ei channels the NNLO is worse than the 
NLO approximation. We emphasize that these curves are 
not obtained from a fit to the data. 



V. SHORT DISTANCE BEHAVIOUR FOR OPE: 
FIXED POINTS AND LIMIT CYCLES 

In this section we analyze the short distance behaviour 
of the equations for the scattering lengths for the singlet 
^5*0, Eq. l(5T|l . and the triplet ^5*1, Eq. ifK^ . channels 
in the short distance limit. According to Eq. H14|l this 
is equivalent to study the mixed boundary condition at 
short distances. 

We study first the case of OPE in the singlet ^5*0 chan- 
nel. At short distances R << the OPE potential 
behaves like the Coulomb potential. Eq. H51() can be eas- 
ily solved in two extreme cases, ao << i? and ao » R. 
While in the first case we get 



ao{R) 



, ttQ « R (72) 



32^/2 

in the second case one solution behaves as 

aoiR) ^ g-l^iM^ iog(fl/flo) ' (^o«R (73) 
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FIG. 1: Evolution of the ^^i, ^Di and Ei NN-threshold parameters from the trivial values at the origin using the Nijmll 
potential to infinity. Top panel: scattering lengths q:o(^) (in fm),ao2(i?) (in fm^) and a2{R) (in fm^). Bottom panel: effective 
ranges vo{R) (in im),vo2{R) (in fm'^) and V2{R) (in fm^). In the left panel we represent a global picture from the Ifm to 40 
fm. The lower panel is a detailed picture in the short distance region below Ifm. 



where Ro is a reference scale fulfilling R < Rq << 1 / rn^ . 
As we see, ao{R) goes to zero in both case but, while 
Eq. (|72) goes rapidly a^iR) 0, Eq. |(72|) goes very 
slowly and with ao(i?) — oo at short distances. In 



momentum space the i? — > limit corresponds to the ul- 
traviolet limit. Eq. H73|) resembles a sort of asymptotic 
freedom and, hence we have an ultraviolet fixed point. 
One can see that the first case, Eq. H72|l . corresponds 



10 




I 1 1 1 1 1 1 1 -J -40 ' ' 1 ' — I ' ' 1 ' 1 

50 100 150 200 250 300 350 400 50 100 150 200 250 300 350 400 



p ( MeV) p ( MeV) 

15 



10 



m 

0! 
W 

r 
Q. 







-5 

50 100 150 200 250 300 350 400 

p(MeV) 



ER-LO 
ER-NLO 
ER-NNLO 




FIG. 2: Theory without exphcit pions. ^^i, ^Di and Ei at LO (contact terms), NLO {k^ terms) and NNLO ( fc"* terms) 
predicted phase shifts for the triplet channel within a pure effective range expansion approximation, with the low energy 
threshold parameters obtained from solving the evolution equations for the threshold parameters, Eq. 14611 with the Nijm II 
potential and using the reduced K-matrix, Eq. 13511 . Data are the PWA from Ref. [32|. 



to selecting the regular solution at the origin, whereas 
Eq. (|75|l corresponds to a generic case, which always con- 
tains an admixture of the irregular solution. The regular 
case at the origin corresponds to integrate from the origin 
starting with the trivial initial condition SJk, 0) = up 
to infinity. As we have discussed in Ref. |3J| the result 
corresponds to a pure OPE interaction, with no short- 
distance interactions. The important thing to realize is 
that regardless of the value of ao at infinity, removing 
one-pion exchange goes into the same value at the ori- 
gin, as implied by Eq. (f7S)l . This also implies that any 



tiny deviation of the ao (R) at small distances results in 
huge variations at infinity. Thus, removing OPE in the 
^5*0 channel implies an extreme fine tuning of the scatter- 
ing length at short distances, and hence of the boundary 
condition at the origin. 

We turn now to the case of the '^Si — Di channel, 
where the tensor force plays a role. In the region close to 
the origin the wave function oscillates wildly and hence 
a WKB approximation may be used. The calculation 
is simplified by taking into account that for the OPE 
interaction the potential matrix is diagonalized by an r- 
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independent unitary transformation, i.e. 



MU(r)M 



ith 



-i_f Uc{r)-WT{r) 







Ucir) + 2UTir) I ' 



V2 1 



(74) 
(75) 



Note that this transformation does not diagonaUze the 
fuU potential U + P/r^ including the centrifugal barrier, 
which for r —f may be neglected. Thus, in the short 
distance limit we may decouple all our equations into 
pairs, and in particular we can apply the transformation 
to the boundary condition, Eq. ||HJ| at zero energy 



MLo(i?)M-i - diagih{R),h{R)) , 



(76) 



where li{R) and hiR) are the logarithmic derivatives 
at zero energy of the decoupled problem with potentials 
Ui = Uc - 4C/t and U2 = Uc + ^Ut respectively. After 
straightforward algebra we get 



MR) 
0:02 (-R) 
a2(i?) 



Rl2iR){Rli{R) + I) - 2 



/i(i?)(3i?/2(i?) + 2) ' 



3 Al2{R) + li{R){iRl2{R) 
R^ h{R){Rl2{R) - I) - 2l2{R) 
15 AhiR) + li{R){mi2{R) + 2) 



(77) 



(78) 



(79) 



Now, as we approach the origin the tensor potential dom- 
inates, and the potential Ui and U2 behave as repulsive 
and attractive Ijr^ potentials respectively, correspond- 
ing to take /i — > cxD and l2{R) by the zero energy limit of 
the logarithmic derivative of a WKB function. 



ao(i?) 

"02 (-R) 

a2{R) 



R 



3RIwkb{R) 
3RIwkb{R) + 
V2R^ 



1 



3 SRlwKBiR) 
R^ RIwkb{R) — 1 



15 3RI 



WKB 



(R) 



(80) 



(81) 



(82) 



with 



As we see the scattering lengths ao, ao2 and a2 present 
on oscillatory behaviour as we approach the origin, so 
they do not converge to a well defined value; as we ap- 
proach to the origin the a's take all possible values. This 
situation corresponds to a limit cycle at short distances. 
A way of avoiding the unbound variation of the scattering 
lengths consists of going to the origin stepwise through 
some envelope subsequence defined by a fixed condition 
for /wkb(-R)- For instance, if we define a cycle by the 
condition a(i?„) = 0, we have RnlwKB{Rn) — 0, yield- 
ing 



(83) 



Here A is an energy independent phase, and Rq a refer- 
ence point, given by 



RolwKBiRo) = l + y^cot{A) , (84) 



and 



Rm = TT?^ = 16fm . 



(85) 



aoiRn) 

a02(-Rn) 
MRn) 



^ 



V2Rl 
6 

30 



Another possibility would to take ZwKB(^ri 
which case one has 



ao{Rn) 

a02(-Rri) 

a2(i?„) 



Rn 



R^ 
30 



As we see, there are infinitely many such possibilities, al- 
though all of them go towards the trivial values, q;o(0+) = 
0^02 (0'*') = 0:2(0+) = 0. Actually, any of the choices cor- 
responds to a different starting condition at infinity, mod- 
ulo a cycle. Conversely, if we go to very short distances, 
where the scattering lengths vary wildly; any tiny per- 
turbation there results in a completely different value at 
infinity. So, we see again that an extreme fine tunning of 
the threshold parameters at short distances is required. 
Finally, let us mention that close to the origin the se- 
quence of cycles can be determined by the solution of the 
equation 



- + XnCOi{/\ + Xn) = 



fl (86) 



in the limit Xn 0. 

In practical numerical calculations the finite integra- 
tion step AR provides a given resolution scale, and these 
infinite limit cycles may not observed due to the rapid os- 
cillations. Instead, one sees the envelope corresponding 
to the stationary points of the scattering lengths. This 
point will become clear below. Sect. IVll 



VI. NON-PERTURBATIVE SOLUTIONS 

A. Evolution of the Low energy parameters 

The exact mathematical analysis of the general set of 
Equations is rather complicated since we are dealing with 
a non- linear system of equations. In Ref. |38j | simple cases 
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are analyzed analytically and the general features which 
can be deduced there are consistent with the numerical 
results we have obtained in the present work. 

As we have said the set of equations, Eq. H51|) and 
Eq. (|58() can be numerically solved. Given the fact that as 
we approach the origin the tensor part of the potential de- 
velops a singularity it is important to carefully check for 
numerical accuracy at short distances. A crucial property 
which must be fulfilled by any algorithm is that of exact 
reversibility; i.e. evolving upwards or downwards should 
be inverse operations of each other. This is an stringent 
test and, moreover, the only way to make sure that when 
the long range piece of the potential is switched on for 
the K-matrix integration we have consistency with the 
effective range expansion up to the relevant order (see 
also below). We prefer to impose this reversibility ex- 
actly, independently on the number of mesh points used 
in the integration, so that any numerical irreversibility 
is merely attributable to computer arithmetic round-off 
errors. This feature will prove extremely relevant when 
computing the phase shifts below since our calculation 
requires upwards integration from lower distances. In all 
calculations presented in this paper we have checked that 
the correct threshold behaviour is obtained. 

Quite generally, we find stable results when we take 
the long distance cut-off to be i?oo = 20fm. On the 
other hand the lowest radius we can achieve numerically 
and preserving reversibility is Rs = O.lfm, mainly due 
to computer arithmetic round-off errors triggered by the 
singularity of the potential. One could further lower the 
radius by a semi-classical approximation as outlined in 
Sect. fVl since as the origin is approached the wave func- 
tion undergoes an increasing number of oscillations and 
WKB methods can be applied. Nevertheless, as we will 
see below, for our short distance cut-off the phase shifts 
for CM momenta up to fc = 250MeV are rather stable 
numerically. 

The strong dependence of the low energy threshold pa- 
rameters on the short distance cut-off provides a clue to 
the fact that there seems to be a lower finite limit for 
the boundary radius Rs = i?min ^ 1.4fm 's^ with still 
an acceptable fit; if the boundary radius is lowered, the 
parameters encoding the short distance boundary con- 
dition which are used as fitting parameters depend in a 
non-smooth way on Rs ■ In addition, the strong singular- 
ity at the origin triggers a fine tuning in those parame- 
ters. According to our previous discussion, this short dis- 
tance fine tuning of low energy parameters is absolutely 
necessary to comply with the independence of the scat- 
tering amplitude on the short distance boundary radius. 
For such a situation, a fit based on successive adiabatic 
changes of Rs becomes impractical since the fitting pa- 
rameters do not change adiabatically and also because 
these parameters should have to be determined to ex- 
traordinary high precision. In addition, the way how the 
limit Rs should be taken differs from channel to 
channel. Our method provides a practical way to over- 
come the difficulty, given the fact that the boundary ra- 



dius is taken exactly to zero along the renormalization 
trajectories while keeping the low energy threshold pa- 
rameters at fixed values. 

1. ^So and Si -without mixing channels 

In Fig. Q we show our results for the evolution of 
the threshold parameters ao, tq and V2 in the singlet 
and triplet '^S'l-without mixing (i.e. neglecting the tensor 
force) channels. The main difference one can appreciate 
from the comparison of both channels is that while the 
scattering length for the ^5*0 channel exhibits a mono- 
tonic trend towards the origin, the scattering length in 
the ^Si channel diverges at a distance of about 0.7fm. 
The interpretation of this fact in our framework is clear; 
the central part of the OPE potential is purely attrac- 
tive. Thus, by eliminating the pions down to a certain 
distance, we are effectively building some repulsion, un- 
til we lose a bound state. An alternative interpretation 
is that as we switch on the OPE potential from the ori- 
gin up to certain distance we can accommodate a bound 
state above 0.7fm. With this interpretation in mind, we 
should add 180*^ to the '^S'l phase shift to comply with 
Levinson's theorem. 



2. '^Si Di channel 

We finally analyze the triplet '^Si — ■^13 1 -channel taking 
into account the tensor mixing. In Fig. Q we show our 
numerical solutions of the set of Eqs. (|46|l . Starting at 
sufficiently long distances (in practice i?oo = 20fm turns 
out to be adequate) and evolve downwards to the ori- 
gin. Operationally this corresponds to eliminate OPE in 
the triplet channel. As can be clearly seen for distances 
above R ~ 3fm nothing dramatic happens and a mono- 
tonic trend is observed. At smaller distances ~ 2fm , 
however, we note a rapid change in the running scatter- 
ing lengths. Again, a rather flat evolution follows until 
the region below Ifm. A zoomed picture is plotted in 
the lower panel for distances shorter than Ifm where the 
cyclic structure of evolution becomes evident, as expected 
from our analysis in Sect. (|V|l . The number of cycles in- 
creases without any bound as the origin is approached. 
This situation is dramatically different from that found 
in the case without tensor mixing, since there OPE pro- 
duced an ultraviolet fixed point. The situation we en- 
counter here is not new and has already been described 
in the context of non-coupled channels. The limit cycle 
structure naturally raises the problem of undefined val- 
ues of the short distance parameters as we take the limit 
R ^ 0. The point is that there is a way of taking the 
limit through equivalent points defined by the property 
aiRn) = ct{Rn+i)', any two such points produce iden- 
tical low energy parameters at infinity. Thus, the limit 
Rn oo through equivalent points produces the same 
parameters at long distances. The cycles in ao2 and a2 
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are hardly seen in the plot due to a low resolution AR 
compared with the typical cycle spacing. 

B. Phase-Shifts 

The standard way of proceeding would be to determine 
the low energy constants or, equivalently, the short dis- 
tance parameters directly from a fit to the data in a large 
energy range ( say up to /c ~ TOtt where the two pion ex- 
change left cut opens up should start contributing) for 
the theory with OPE. The low energy parameters would 
have to be recomputed, and the description at lower en- 
ergies {k < rriT, ) might become even worse than a pure 
effective range expansion (see e.g. Refs. 39, 40]). Obvi- 
ously, this is an undesirable situation. The effective range 
expansion is convergent up to the OPE left cut, located 
at fc — ±im^/2 and should be applied only there ^. Our 
formalism can be specifically constructed to avoid such a 
situation. Once the threshold parameters are determined 
in the short distance limit Rs 0, our phase shifts be- 
come pure predictions without any additional parameter 
fitting obtained to a given order fc^ expansion of the ini- 
tial condition by integrating Eq. 148|l using the effective 
range type of initial condition, 

Ks = K(i?5) = -a^i + ^rsk^ + vgk^ + ... (87) 

with with i?s ^ 0. The solution of Eq. (gHJ at R 
oo gives a solution which when expanded in powers of 
k^ exactly reproduces ERE to the order imposed by the 
initial condition, Eq. H35() . Thus, the difference beyond 
the displayed terms is merely attributable to the OPE 
potential. 

In what follows we use LO, NLO,NNLO, etc. to de- 
note keeping up to the first, second, third order terms in 
Eq. l|S7jl respectively. 

1. ^ So and ^ Si-without mixing channels 

In Fig. lO we show the results for the phase shifts 
for both ^5*0 and ^5*1 -without mixing channels depend- 
ing on the number of terms kept in the low energy ex- 
pansion at short distances. Our results exhibit a good 
convergence rate. For comparison we also depict the 
effective range expansion results without explicit pions, 
which is expected to work at low energies only. As we 



* The fact that only two terms in the expansion, involving the 
scattering length and the effective range only, works so well 
at high momenta, almost up to fc ~ 77i,r, is purely accidental. 
Actually, including the next V2 term in the expansion and fit- 
ting it in the region k < r7i^/2 does not reproduce the data for 
m,r/2 < k < m^, but improves the fit for k < m,r/2. This is ob- 
viously an indication of the breakdown of the expansion beyond 
the analyticity domain. 



see, the effect of introducing pions always improves the 
results. This can be fully appreciated at NNLO, where 
ERE does a poor job above CM momenta lOOMeV, 
but explicit OPE effects enlarge the energy range up to 
about ^ 140MeV ^ to^. where we expect explicit two 
pion exchange contributions to start playing a role. 

2. ^Si - 3Di channel 

Once the short distance evolution of the low energy pa- 
rameters are known one may compute the phase shifts to 
any order of the approximation in a fc^ expansion of the 
initial condition without any additional parameter fitting 
by integrating Eq. 148(1 upwards with a suitable initial 
condition at a short distance radius. As a matter of fact 
the practical choice of the radius in the numerical cal- 
culation is far from obvious, particularly in the triplet 
channel case where the low energy parameters take un- 
bounded values in an increasingly finer scale at short dis- 
tances (see e.g. Fig. It is most practical to use the 
WKB approximation to match the numerical solution at 
a radius i?wKB which can safely be taken in the range 
~ 0.5fm. The results for LO (contact terms), NLO (fc^ 
terms) and NNLO ( k^ terms) are presented in Fig. B 
and compared to the partial wave analysis of Ref. |32 |. 
As we see the best scheme to take into account the OPE 
potential corresponds to use the NLO initial condition. 
This means on the one hand that while the scattering 
lengths may be considered large and comparable to the 
effective ranges the curvature parameters V2 can be con- 
sidered to be small. 



C. Finite cut-off effects 

Finite short distance cut-off effects in the scattering 
phase shifts can be seen in Fig. |H1 for finite radii Rs = 
1.4fm and Rs = l.Sfm as compared to the renormalized 
i?5 = case, for the OPE-LO, OPE-NLO and OPE- 
NNLO approximations. As one naively expects these fi- 
nite effects increase for larger energies, since they probe 
smaller wavelengths. A very important feature which can 
be deduced from the plots is that these effects are sizable 
for momenta where TPE effects should not play a deci- 
sive role TOtt/Z < fc < m,r role. Thus, letting a finite 
short distance boundary radius Rs ''^ 1.4fm provides a 
large systematic error, already in the region where OPE 
dominates. Thus, it is not clear whether TPE can bee 
seen in the central NN waves with a finite cut-off distance 
of about Rc — 1.4fm. Of course, one should include TPE 
contributions in order to make a definite statement. 



D. Are pions perturbative ? 

The discussion of which power counting is the appro- 
priate one for the NN interaction corresponds physically 
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FIG. 3: Evolution of the scattering (left panel) and ^^i-without mixing (right panel ) NN-threshold parameters ao{R) (in 
fm),ro(i?) (in fm) and V2{R) (in fm^) from the asymptotic values at infinity (which we take in practice Roo = 20fm ) when 
OPE effects are removed down to the origin. 



to the question whether or not the pion cloud can be con- 
sidered to be perturbative. It is important to realize that 
within our framework we are considering OPE depar- 
tures from the effective range expansion to a given order. 
Thus, at sufficiently low k explicit pion effects can al- 
ways be considered perturbative. This is so regardless of 
the number of terms included in the initial condition. 
Actually, the point is rather if the low energy threshold 
parameters can be considered large or small. According 
to our results in Fig. © it seems that best possible agree- 
ment can be obtained when both the scattering lengths 
and the effective ranges are taken to be large, while other 
low energy parameters can be taken to be small. To prop- 
erly emphasize this point we plot in Fig. Q the scaled 
K matrix computed including OPE and compared to the 
ERE to LO, NLO and NNLO. Given this fact we expect 
a kind of consistent long distance perturbation theory to 
work. The details of such an expansion will be presented 
elsewhere. 



E. Evolution of the short distance boundary 
condition 

As we have said, the short distance singularity of the 
OPE potential enforces a very precise determination of 
the running low energy threshold parameters at short 
distances, and hence of the boundary condition. We can 
directly determine this dependence by using Eq. H62|l and 
Eq. Ht)5|) . For simplicity and to illustrate the point we 
just display in Fig. (PJ the behaviour of the boundary 
condition parameters as a function of the short distances 



boundary radius in the zero energy limit, both for the sin- 
glet ^Sq and triplet ^5*1 channel without mixing and for 
the triplet — ^ Di channel. The fixed point and limit 
cycle behaviour obtained for the running of the low en- 
ergy threshold parameters maps into a similar behaviour 
for the short distance boundary condition. From the pic- 
ture it is clear that the standard procedure of integrating 
the Schrodinger equation upwards from a given short dis- 
tance boundary radius to infinity in order to fit the low 
energy parameters would require a very high precision 
determination of a rapidly varying boundary condition 
in the case of the triplet ^5*1 — ^ Di channel. It is clear 
that a determination of the Cq constants from a fit to the 
phase shifts in the low energy region would be extremely 
delicate in the limit Rc —^ in practice. Instead, the 
present approach computes directly the boundary con- 
dition in a power expansion of the energy at any given 
radius from the physical values of the low energy parame- 
ters. Actually, our method is equivalent to integrate the 
Schrodinger equation from that short distance bound- 
ary radius to infinity. In addtiiton, the singular and 
attractive nature of the OPE potential allows a WKB 
treatment of the short distance singularity, and allows to 
eliminate the finite cut-off radius taking the limit i?c — > 0. 
Obviously, the present framework can be extended to re- 
analyze the role of TPE potentials in a non perturbative 
way and completely free of finite cut-off artifacts, where 



15 




FIG. 4: Evolution of the "^Si, '^Di and Ei NN-thrcshold parameters from the physical values at infinity down to the origin 
using the OPE potential. Top panel: scattering lengths Qo(-R) (in fm),ao2(i?) (in fm^) and a2{R) (in fm^). Bottom panel: 
effective ranges vo{R) (in fm),vo2{R) (in fm^) and V2{R) (in fm^). In the left panel we represent a global picture from the Ifm 
to 20 fm. The lower panel is a detailed picture in the short distance region below Ifm. Limit cycles are clearly visible in the 
s-wave scattering length ao and effective ranges. The Ei and ^Di scattering lengths ao2 and a2 go quickly to zero below 0.25 
fm. 
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FIG. 5: Predicted phase shifts in the ^So (left panel) and the ^S'l-without mixing (right panel) channels for NN scattering as 
a function of the CM momentum in MeV. In the '^Si channel we assume no mixing according to Eq. H48^ when OPE potential 
is switched on and the initial condition is a low energy expansion of the A'-matrix at short distances, (See Eq. 18711 in the 
main text). LO means keeping ctsfi only, NLO keeping as,o and ro,s and NNLO keeping as,o, '"o.s and V2,s- The short range 
parameters are directly determined by evolving the low energy parameters from their experimental values, Eg. 166^ and Ea. l67|l 
ER-LO, ER-NLO and ER-NNLO corresponds to a pure effective ran ge e xpansion keeping ao only, as,o and ro, ao, '"o and V2 
respectively. No further fit is involved. Data are the PWA from Ref. |32 |. 



the short distance behaviour is qualitatively similar ^ 



VII. CONCLUSIONS 

In the present paper we have analyzed the renormal- 
ization of the OPE interaction in the presence of contact 
and derivative interactions of any order for NN scatter- 
ing both for the singlet and triplet channel states. The 
basic point of our approach is to regularize the unknown 
short distance physics by means of a boundary condition 
at a certain boundary radius, above which the OPE po- 
tential is assumed to work, i.e. where pions are treated 
explicitly. Below that scale pions contribute implicitly to 
the scattering properties although always in combination 
with other effects which cannot be disentangled unless 
a given distance scale is specified. Actually, when the 
boundary radius goes to infinity, above the pion Comp- 
ton wavelength, we have a low energy theory of contact 
interactions and derivatives there-off. As the boundary 
radius goes below the OPE range, we have a theory where 
pions are eliminated above the scale set by the boundary. 



^ Unlike the OPE where one has both an attractive and repulsive 
singularity (see Sect.EJ, in the TPE case one encounters 
attractive singularities for coupled channels. 



This allows to remove explicitly pion effects in the thresh- 
old parameters for the OPE potential in an unambiguous 
and model independent way. The renormalization group 
flow implied by our non-perturbative equations is unique 
provided the OPE potential is assumed to be valid al 
the way down to the origin. This is obviously not a re- 
alistic assumption but it is absolutely necessary to go 
to these small distances in order to get rid of any finite 
short-distance cut-off effect and properly define the OPE 
contributions to the scattering observables. This result 
fully complies to the spirit of an effective field theory, 
the terms in a low momentum expansion of the ampli- 
tude are shape independent while the remaining powers 
depend both on the long distance OPE details (like the 
left branch cut) and the shape independent low energy 
parameters themselves. 

The short distance behaviour of threshold parameters 
present either a ultraviolet fixed point structure in the 
^5*0 and 5*1 -without-mixing channels whereas we find 
limit cycles for the '^S'l Di channel due to the singular 
and attractive behaviour of the OPE contribution to the 
tensor potential. This means that in the latter case there 
is not a monotonous trend at short distances. A direct 
consequence of having both ultraviolet fixed points and 
limit cycles for the threshold parameters is that a deli- 
cate fine tunning of the short distance physics is implied. 
In addition, for the experimental values of the threshold 
parameters one obtains huge changes for distances below 
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FIG. 6: ^Si, and E\ at LO (contact terms), NLO {k^ terms) and NNLO ( k"^ terms) predicted phase shifts for the triplet 
channel. The initial condition is a low energy expansion of the if-matrix at short distances. By construction the low energy 
parameters a , r and v coincide with those extracted from the Nijmll potential. Data are the PWA from Ref. |32j. 



2fm when OPE effected are removed. Nevertheless, we 
find moderate changes in the phase shifts due to expUcit 
pion effects. Actually, in the ^Sq and ^S'l -without mix- 
ing channels the effect is found to be compatible with 
a perturbative treatment. In the '^S'l Di channel the 
effect is a bit more complicated due to the presence of ul- 
traviolet limit cycles triggered by the singular character 
of the tensor potential; the coupled channel amplitudes 
are non-perturbatively renornializable while they become 
perturbatively non-renormalizable. This makes a naive 
perturbative treatment slightly more subtle. One of the 
advantages of having a renormalizable theory is that non- 



perturbative equations make sense, and any perturbative 
treatment should arise as a controllable approximation to 
the full equations. As we have pointed out along the pa- 
per, this is probably an advantage of using coordinate 
space methods and a boundary condition renormaliza- 
tion versus momentum space methods. 

Taking into account all the nice features of the present 
calculation, in particular, getting a handle on the finite 
cut-off corrections, the results presented in this paper are 
very satisfactory suggesting several improvements. Ex- 
plicit Two Pion Exchange contributions are expected to 
contribute significantly at about 1.5 — 2fm at the level 
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FIG. 7: The effect of including the OPE potential on top of the effective range expansion (ERE). Top panel : '^Si, '^Di and 
El at LO with (OPE-LO) and without (ER-LO) OPE explicit effects. Middle panel: same but for NLO. Bottom panel: same 
but NNLO. The difference between ER and OPE indicates the size of the explicit effects due to the OPE potential to LO 
(contact terms), NLO (k^ terms) and NNLO ( fc* terms). In both cases the low energy threshold parameters coincide with 
those extracted from the Nijmll potential. Data are the PWA from Ref. [s^ . 



of the potential, so our results for the evolution of the 
threshold parameters should not be considered realistic 
below that scale, or equivalently above CM momenta of 
about 100 — ISOMeV, as it seems to be the case. In addi- 
tion, our description should be enlarged to include higher 
partial waves. For peripheral waves one expects pertur- 
bative methods to work since there is a strong centrifugal 



suppression of the wave function at the origin, and per- 
turbative renormalization methods can be applied. For 
those the present approach does not have much to say. 
Low partial waves, however, are particularly interesting 
since a re-summation of pion exchanges seems crucial to 
understand the data. Work along these lines will be pre- 
sented elsewhere [38|. 
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FIG. 8: The effect of having a finite short distance boundary radius for the OPE potential on top of the effective range expansion 
(ERE). We compare the theory with finite radii R = 1.4fm and R = 1.8fm with the renormalized theory R = 0. Top panel 
: ^5*1, ^-Di and Ei OPE-LO. Middle panel: same but for NLO. Bottom panel: same but NNLO. In all cases the low energy 
threshold parameters coincide with those extracted from the Nijmll potential. Data are the PWA from Ref. [s^- 

Acknowledgments APPENDIX A: ANOTHER DERIVATION OF 

THE VARIABLE S-MATRIX 



We thank J. Nieves for discussions. This work is sup- 
ported in part by funds provided by the Spanish DGI 
with grant no. BMF2002-03218, Junta de Andalucfa 
grant no. FM-225 and EURIDICE grant number HPRN- 
CT-2003-00311. 



In order to deduce a variable S-matrix equation, we 
determine first the infinitesimal change of the 5' matrix 
under a general deformation of the potential U(r) — s- 
U(r) + AU(r). Using Schrodinger's equation (O and the 
standard Lagrange's identity adapted to this particular 



20 





FIG. 9: Evolution of the dimensionless short distance boundary conditions at zero energy Co(-R) = 1 — -RLo(-R) = 1 — 
Ru'q{R)uo{R)~^ with the boundary radius R due to OPE potential. Top: Singlet ^So channel (left) and Triplet ^Si channel 
(right) without mixing using Eq. ^ . Bottom: The triplet ^Si Di channel. The coefficients and C^^ are related 

to the running scattering lengths aoo , Q.02 and 022 through Eq. (I65II . Large scale (left) and zoomed (right) picture. 



case, we get after integration In particular, for the parametric family of potentials 

U(r, R) =. e{R - r)U(r) we get 

[u(r)tAu'(r)-u'(r)tAu(r)] |J = y^ dm(r)t AU(r)u(r) 2ikS\R)S' [R] = n{R)^\J{R)n{R) (A3) 

(Aland using the value of the wave function at the outer 
which, after using the asymptotic form of the matrix wave boundary 

function, Eq. ®, yields u(i?) = \i^-\R) - h^+^R)S{R) (A4) 

2ifcS^AS - r dru(r)^AV(r)u(r) (A2) ^""^^^^ ^""^ ^^"^ variable S-matrix equation, Eq. CEl. 

^ ' \ J \ J \ J Note, that the variation of the potential is done with a 
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fixed boundary condition. 



APPENDIX B: DIFFICULTIES IN EXTRACTING 
THE LOW ENERGY PARAMETERS 

In this appendix we want to elaborate on the prob- 
lems we have encountered while fitting the NN data base 
of Ref. [s^ within a generalized coupled channel effective 
range expansion, Eq. f l35|l . Unfortunately, this data base 
does not provide error estimates for their phase shifts (al- 
though 8 significant digits), nor the typical energy reso- 
lution where these data should be trusted, so some com- 
promise must be made. 

We use the NN-data and define the as 

where we take Ai?LAB = O.OlMeV , and .R'nn and AXnn 
are the mean value and the standard deviation of the six 
potentials listed in the NN-data base [22|j which can be 
taken as independent uncorrelated primary data. The 
factor Mj (4p) is the Jacobian of the transformation be- 
tween the Lab-energy and the CM. momentum, i?LAB = 
2p^/M, and would correspond to make an equidistant 
sampling in p, in the limit AE'lab (this is why we 
take a small energy spacing) . This weight factor is intro- 
duced in order to enhance the region at low momenta. 
On the other hand, very low momenta, must be excluded 
since the resulting mean value _fC-matrix is incompatible 
within the attributed errors with the expected theoretical 



behaviour, Eq. (|2SJ, thus we take i?LAB > O.SMeV Also, 
the fit goes up to £^lab < lOMeV, which corresponds to 
a CM. momentum about p = 'rn-^/2 where we expect the 
finite polynomial of the scaled K — matrix to truly repre- 
sent an analytical function within the convergence radius 
up to the branch cut singularity located aX p — ±im7r/2. 
The form of the fitting function is 

-ft'ER = -/9 + ro/ -I- W2/ + W3/ + W4/ + . . . (B2) 

In Fig. IIUI we show as an illustration the V2 parameter 
determined from a fit to the low energy region of the NN 
data base as a function of the maximal LAB-energy 
considered in the fit. As we see, instead of a plateau 
within some energy window, we observe an ever changing 
value. We observe no stability depending on the number 
of terms considered in Eq. ljB2|l either. For comparison 
we also plot the values we obtained by integrating the 
Eqs. l|nU,(|n2l and with the Nijmll and Reid93 po- 
tentials in Sect. lIVI which where quite stable numerically. 
As we see, the values obtained from the fit, in the cho- 
sen energy window are hardly compatible. The deceptive 
features extend to other channels, and non diagonal low 
energy threshold parameters such as the matrix elements 
of a and r. 

Finally, we have also tried, with no success, other 
methods for the determination of the low energy thresh- 
old parameters, like evaluation of derivatives within sev- 
eral algorithms. The reason for the failure has to do with 
round-off errors generated by the relatively small number 
of digits provided in the NN database. 
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